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Radiation Mechanism of the Soft Gamma-ray Pulsar PSR B 1509-58 
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The outer gap model is used here to explain the spectrum and the energy dependent light curves 
of the X-ray and soft 7-ray radiations of the spin-down powered pulsar PSR B1509-58. In the 
outer gap model, most pairs inside the gap are created around the null charge surface and the 
gap's electric field separates the two charges to move in opposite directions. Consequently, the 
region from the null charge surface to the light cylinder is dominated by the outflow of particles 
and that from the null charge surface to the star is dominated by the inflow of particles. The inflow 
and outflow of particles move along the magnetic field lines and emit curvature photons, and the 
incoming curvature photons are converted to pairs by the strong magnetic field of the star. These 
pairs emit synchrotron photons. We suggest that the X-rays and soft 7-rays of PSR B 1509-58 result 
from the synchrotron radiation of these pairs, and the viewing angle of PSR B 1509-58 only receives 
the inflow radiation. The magnetic pair creation requires a large pitch angle, which makes the 
pulse profile of the synchrotron radiation distinct from that of the curvature radiation. We carefully 
trace the pulse profiles of the synchrotron radiation with different pitch angles. We find that the 
differences between the light curves of different energy bands are due to the different pitch angles 
of the secondary pairs, and the second peak appearing at E > lOMeV comes from the region near 
the star, where the stronger magnetic field allows the pair creation to happen with a smaller pitch 
angle. 



I. INTRODUCTION 

PSR B1509-58 (hereafter PSR B1509) is a unique 
subject of the multi-wavelength study of the high en- 
ergy radiation mechanism of the spin-down powered 
pulsar. It is a young (about 1600 years old) spin-down 
powered pulsar with a period of about 150 ms and a 
surface magnetic field of ~ 1.5 x 10 13 G, which is much 
stronger than those of most of the canonical pulsars 
(~ 10 12 G). Multi-wavelength observations show that, 
from 0.1 keV to 10 MeV, PSR B1509 has a wide sin- 
gle peak dHl], whose shape can be described by two 
Gaussian functions, when E > lOMeV, anotherpeak 
shows up near the peak of radio emission 0, S • 

PSR B1509 is unique for its soft 7-ray spectrum and 
a very strong surface magnetic field among the known 
7-ray spin-down powered pulsars. It is suggested that 
there may be a new class of spin-down powered pul- 
sar with high magnetic field (>10 13 G) and soft 7-ray 
radiation, called soft 7-ray pulsar |8f. And there is 
at least another pulsar that is similar to PSR B1509, 
PSR J1846-0258, a young pulsar with P ~ 324ms and 
surface magnetic field B s ~ 4.9 x 10 13 G, which cannot 
be detected by Fermi but by X-ray telescopes @, [Io| • 
Therefore, a multi-wavelength study of the radiation 
mechanism of the PSR B1509 can be a bridgehead to 
investigate this possible new class of pulsars. 

Recent observations show that, two kinds of pulsars: 
the spin-down powered pulsar with B s ~ 10 12 G and 
the magnetar with B s ~ 10 14 — 10 15 G, have no clear 
boundary between them. It is believed that the radi- 
ation of a magnetar (AXP or SGR) is powered by the 
energy stored in the super strong magnetic field with 
B s > 10 14 G [Tll| - However, the discoveries of the low 
magnetic field soft 7-ray repeater [T2I . [iH ] , the radio 
emission from magnetar 61 ] and the magnetar-like 



outburst of spin-down powered pulsar with high mag- 
netic field |9j challenge the existing theoretical mod- 
els. The study of the radiation mechanism of the very 
high magnetic field spin-down powered pulsar, and the 
connection between the canonical spin-down powered 
pulsar and magnetar, can help us to understand the 
physics of the two groups of neutron stars. As a family 
member of the high magnetic field spin-down powered 
pulsars, PSR B1509 may provide some insights for un- 
derstanding the differences between typical spin-down 
powered pulsars and magnetars. 

Here we give explanations for the spectrum and 
the energy dependent light curves of the X-ray and 
soft 7-ray radiations of the spin-down powered pul- 
sar PSR B1509-58, based on the outer gap model 
(l7| . Cheng, Ruderman & Zhang 18] showed that 
most electron/positron pairs created inside the outer 
gap are produced near the null charge surface; the 
strong electric field inside the gap separates the op- 
posite charges to move in opposite directions. There- 
fore from the null charge surface to the light cylinder 
the radiation is mainly outward, and from the null 
charge surface to the star the radiation is mainly in- 
ward. Since the light cylinder of PSR B1509 is much 
bigger than that of the Crab pulsar, most outflow cur- 
vature photons can escape from photon-photon pair 
creation. If the line of sight is in the direction of out- 
going radiation beam, the spectrum of PSR B1509 
should be a characteristic pulsar spectrum, namely a 
power law with exponential cut-off spectrum as pre- 
dicted by Wang, Takata & Cheng [Ti| . Obviously this 
conflicts with the observation. Therefore, we propose 
that the viewing angle of PSR B1509 is in the direc- 
tion of incoming radiation beam, where the curvature 
photons of hundreds MeV are converted into pairs by 
the strong magnetic field. The observed X-rays and 
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FIG. 1: The sketch of the structure of the gap. 



soft 7-rays are synchrotron radiation of these pairs. 

II. SIMULATIONS OF THE SPECTRUM 
AND THE ENERGY DEPENDENT LIGHT 
CURVES 

Figure [T] shows the structure of the acceleration re- 
gion. Outside the null charge surface, the structure 
of t he g ap can be simplified as a two-layer structure 
[l9l [20j . From the null charge surface to the stellar 
surface, the accelerating electric field is screened out. 
Outside the null charge surface, the particle can get 
the equilibrium state that the particle loses the energy 
obtained from the potential drop via curvature radi- 
ation. Hence its Lorentz factor can be obtained from 
eE\\{r)c = 2e 2 cy 4 (f)/3s 2 (r). On the other hand, for 
the particle moving from the null charge surface to the 
star, the energy loss by the curvature radiation can- 
not be compensated for by the reduced electric field. 
Therefore, we solve numerically the evolution of the 
Lorentz factor, which is given by 

We trace the field lines in the gap from the stellar 
surface to the light cylinder. At each step f in the 
field line, we calculate the pulse phase and viewing 
angle (if}, £) of the curvature radiation, whose direc- 
tion is tengential to the field line in the co-rotating 
frame. Then we trace the direction of each incoming 
curvature radiation to find the place where the mag- 
netic pair creation happens. The synchroton radiation 
of the pair is a hollow cone, where the (i/j, C) of each 
direction is calculated. 

If the radiation satisfies |£ — /3| < 0.5°, where /3 is 
our viewing angle, we calculate its spectrum. Because 
some synchrotron photons approach the star and be- 
come pairs, we trace each direction of the synchrotron 
radiation to calculate the attenuation caused by the 
absorption of the magnetic field. The simulated phase 
averaged spectrum is the sum of the visible survival 
synchrotron radiation, which is shown by Figure [2J 
The energy dependent light curves are obtained by 
integrating the phase resolved spectra. The number 
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FIG. 2: The simulated phase averaged spectrum of PSR 
B1509 (solid line), comparing with the observed data pro- 
vided by EGERT (triangles), COMPTEL (circles H), AG- 
ILE (filled squares 0), Fermi (diamonds 0), ASCA (dot 
line [Si), Ginga (dot-dashed line [II]), OSSE (long-dashed 
line 0), Welcome (dot-dot-dashed line [H), and RXTE 
(short-dashed line 0]). 




FIG. 3: The simulated energy dependent light curves in 
the pulse phase and energy plane. The color represents the 
percentage of the number of the photons of certain pulse 
phase interval, in the total number of photons of certain 
energy range. 

of the photons with E\ < E < E2 measured at pulse 
phases between ip\ and ip 2 is calculated from 

N^{E X ,E 2 ,^ 2 ) oc f 2 F tot (£,Vi<V <^2)dE. 

J E x 

(2) 

By fitting the observed spectrum and the energy de- 
pendent light curves, we choose the inclination an- 
gle as a = 20°, and the viewing angle as (3 = 11°. 
Figure [3] is the simulated energy dependent light 
curves, and the comparisons between the simulated 
light curves and the observed ones are given by Fig- 
ure m 
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FIG. 4: The simulated energy dependent light curves (grey histograms), comparing with the observed data (solid lines), 
which are provided by Fermi @|, AGILE [ff, COMPTEL [H, OSSE [| and BeppoSAX 



III. DISCUSSION 

• Where are the GeV curvature photons? 

The magnetosphere of PSR B1509 can generate 
GeV 7-ray photons by accelerating the charged 
particles in the outer gap, as well as other 7- 
ray spin-down powered pulsars. However, these 
outgoing curvature photons are missed by our 
line of sight (Figure [SJ. As shown in Figure [SJ if 
the viewing angle or inclination angle increases, 
the GeV curvature photons become visible. 

• Why is the cut-off energy of sepctrum so low? 

This is due to the absorption of the magnetic 
field. The magnetic field of PSR B1509 con- 
verts not only the curvature photons, which are 



emitted by the incoming particles, to pairs, but 
also part of the synchrotron photons emitted by 
these pairs. Figure [7] is the spectra of a = 20° 
and /3 = 11°, where the outgoing curvature ra- 
diation and the possible survival incoming cur- 
vature radiation are missed by the line of sight. 
The dashed and solid lines are the cases with 
and without the attenuation of the synchrotron 
radiation of the secondary pair caused by the 
magnetic field. 

Why is there a second peak at phase~0? 

The second peak of PSR B1509 is the syn- 
chrotron radiation emitted by the pairs with 
smaller pitch angles, which are converted from 
the curvature photons with higher energy emit- 
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FIG. 5: The skymap of the viewing angle £ and the pulse 
phase ip of the outgoing curvature photons (solid line) 
emitted from the null charge surface to the light cylin- 
der, the incoming curvature photons (dashed line) emit- 
ted from the null charge surface to the stellar surface, and 
the synchrotron photons (color) emitted from the places 
where the incoming curvature photons become pairs. The 
inclination angle is a — 20°. 
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FIG. 6: The simulated phase averaged spectra of differ- 
ent inclination angles a and different viewing angles /?. 
The dashed lines are the synchrotron radiations of the 
secondary pairs, the dot-dashed lines are the outgoing cur- 
vature radiations, the dot lines are the survival incoming 
curvature radiations, and the solid lines are the total spec- 
tra. 

ted from the region of 20R S < r < 30R S (Fig- 
ure |SJ). For a fixed birth place of the curva- 
ture photons r, the photons with higher E cur 
become pairs earlier. These photons with higher 
energy can become pairs under smaller pitch an- 
gle P , and the pairs are generated further from 
the star. If a synchrotron photon is emitted fur- 
ther from the star, it has higher chance to sur- 
vive from the absorption of the magnetic field. 
When j3 < a, the pulse phase of the synchrotron 
radiation is closer to when the pitch angle is 
smaller, and closer to 0.5 when the pitch angle 
is larger. 

• Why does the second peak show up when energy 
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FIG. 7: The spectra of the synchrotron radiations of the 
secondary pairs generated by pair creations. The solid 
line is the one without considering the attenuation of the 
synchrotron photons by the magnetic pair creation. The 
dashed line is the spectrum of the survival synchrotron 
photons from the magnetic pair creation. The dot-dashed 
line is the contribution of the photon-photon pair creation, 
which considers the attenuation caused by the magnetic 
pair creation. 




FIG. 8: The energy dependent light curves, whose original 
curvature photons are generated in five regions: r < 20R S , 
20R S < r < 30R S , 30R S < r < A0R s , A0R s < r < 50R S , 
and r > 50R S . 



increases? 

If the particle is made via the magnetic pair cre- 
ation, the typical energy of its synchrotron radi- 
ation is proportional to the energy of its original 
curvature photon. The curvature photon with 
higher energy can become pair under smaller 
pitch angle, which makes the pulse phase of the 
synchrotron radiation close to 0. Therefore, for 
the same birth place of the curvature photons, if 
the viewing angle is smaller than the inclination 
angle, the observed spectrum of the synchrotron 
radiation, whose pulse phase is closer to 0, has 
higher cut-off energy. 

The viewing angle /? determines the energy 
where the second peak shows up. If /3 < a, the 
second peak shows up at higher energy when the 
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P gets closer to a, because the observed second 
peak is contributed by the pairs with smaller 
pitch angles 6 p , which were the curvature pho- 
tons of higher energy. 



IV. CONCLUSION 

The line of sight of PSR B1509-58 is in the di- 
rection of incoming beam instead of outgoing beam, 
otherwise a characteristic power law with exponen- 
tial cut-off spectrum with cut-off energy around a few 
GeV should be observed. In order to avoid seeing the 
outgoing flux and fit the observed multi-wavelength 
light curves we need to choose inclination angle=20° 
and viewing angle=ll°. The observed spectrum is 
the synchrotron radiation emitted by the pairs pro- 
duced by the magnetic field that converts the major 



part of the incoming curvature photons. The energy 
of the curvature photon and the magnetic field deter- 
mine the pitch angle of the pair. We find that the 
differences between the light curves of different en- 
ergy bands are due to the different pitch angles of 
the secondary pairs, and the second peak appearing 
at E>10MeV is the synchrotron radiation emitted by 
the pairs with smaller pitch angles. 
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